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Abstract 

We attempt to describe soft hadron interactions in the framework of saturation models, 
one based upon the Bahtsky-Kovchegov non-hnear equation and another one due to Golec- 
Biernat and Wiisthoff. For pp, Kp, and np scattering the relevant hadronic wave functions 
are formulated, and total, elastic cross-sections, and the forward elastic slope are calculated 
and compared to experimental data. The saturation mechanism leads to reasonable repro- 
duction of the data for the quantities analyzed, except for the forward elastic slope, where 
the predicted increase with energy is too moderate. 
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1 Introduction 

Understanding the high energy behaviour of total hadronic cross-sections within the framework 
of QCD is one of the intriguing problems of high energy physics. The main difficulty lies in the 
fact that presently most applications of QCD are based on perturbation theory which is only 
applicable for "hard" processes (i.e. it needs a "hard" scale), while hadronic processes near the 
forward direction are "soft" and non-perturbative by definition. On the other hand, the past 
few years have seen much activity in the successful apphcation of QCD to DIS processes. For 
values of > 2 GeV^, the use of perturbative QCD seems to be trustworthy. For very small 
one can rely on Regge theory (e.g. PQ) which provides a reasonable description of the data. 
The construction of a very promising bridge between these two theoretical frameworks has been 
pioneered by the concept of high parton densities and saturation. Models based upon this idea 
have been successful in describing the DIS cross-section for all values of and energies x < 



The goal of this paper is an attempt to apply the dipole picture and the physics of high parton 
densities to soft hadronic cross-sections. We want to explore to what extent the high energy 
hadronic asymptotic behaviour can be explained by the saturation hypothesis, which - so far - 
has been tested in the context of deep inelastic scattering at small x and in 7 — 7 scattering. 
The idea of saturation concerns the interactions between partons from different cascades, which 
in the linear evolution equations (DGLAP and BFKL) are not included, and which become 
more important with increasing energy. The parton saturation phenomenon then introduces a 
characteristic momentum scale Qs{x), which is a measure of the density of the saturated gluons. 
It grows rapidly with energy, and it is proportional to ^ 00 EI with A ~ 0.2. Parton saturation 
effects are expected to set in at low values of and x, where the parton densities are sufficiently 
large. 

At this stage we do not discern how to relate the dipole picture to the additive quark model which 
has been successful in explaining and relating different hadronic total cross sections. Instead, 
we consider this study as being exploratory, and we will not attempt to draw any conclusions 
concerning this rather fundamental issue. 

The basis of our endevour is the successful fit |1] to the F2 structure function data for all values 
of and x < 0.01, within the framework of QCD, achieved by using an approximate solution to 
the Balitsky-Kovchegov (BK) |H| nonlinear evolution equation, and adding a correcting function 
to improve the DGLAP behaviour at large Q^. Although soft physics is not explicitly included, 
agreement with experiment is found for all parameters associated with F2, in particular for the 
logarithmic slope A = c}lnF2/9(ln l/x), a value of A ~ 0.08 was obtained at very low x and 
well below 1 GeV^, i.e. in the saturation region. This agrees with the value of the intercept of 
the "soft" Pomeron, associated with the Donnachie-Landshoff (DL) model 9J. 

In this paper we start from the hypothesis that also in hadron-hadron scattering at high energies 
color dipoles might be the correct degrees of freedom, even when large transverse distances come 
into play. We start from the well known expression for DIS cross-sections 
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where denotes the virtual photon's four momentum squared, Vt.l its wave function, W'^ 
the energy squared in the photon-proton system and z, (1-z) the momentum fraction taken by 
the quark (antiquark) respectively, is the transverse distance between the q and g, and 
^ = [w'^+ Q'^) • There are two main elements in Eq. p.l|) : (a) the wave function of the virtual 
photon, and (b) the dipole cross-section, which describes the interaction of the qq with the proton 
target, through the exchange of a gluon ladder. In the DIS case the wave function for the virtual 
photon is well known, whereas for the hadron case this is not so. We discuss the question of the 
hadronic wave functions in Section 3. 

It would be naive for us to expect that our treatment is able to yield the complete hadronic 
cross-section. We have, at least, two reasons for this statement: first, at large impact parameter 
we have to include the non-perturbative contributions even for the so called "hard" processes 
IinilIIllI2] since this behaviour is defined by the spectrum of hadrons [T3|; second, at present the 
impact parameter dependence of the interaction is only treated approximately. In j2] , the dipole 
cross section has a built-in sharp cutoff in h at the value of the proton radius; in other cases j3] 
the equation is first solved for b = 0, and then an ansatz is made regarding factorization and the 
assumed b dependence of the dipole cross section. For hadronic interactions the impact parameter 
dependence is known to be important, and neither a sharp (in particular: energy independent 
cutoff) nor the method of calculating saturation at b = 0, and assuming that the 6-shape does 
not change with energy, is, at best, a very rough approximation to the physical situation. 

The content of the paper is as follows. In Section 2 we discuss the numerical solution of the BK 
equation |8j, and the changes that must be made to adapt this for the calculation of hadron- 
hadron cross-sections. In Section 3 we present the details of our calculation. Section 4 is devoted 
to the overall picture, including comparison of the model predictions with experimental data. 
Section 5 contains a discussion of our results and our conclusions. In the Appendix we explain 
why the results for the forward elastic slope are so shallow. 



In ^ an approximate solution to the BK non-linear evolution equation [S] was obtained using 
numerical techniques. Below we briefly review the method used and the main results obtained. 
For more details of the method of solution we refer to P|. 

The solution of the BK equation which we denote by iV, takes into account the collective phenom- 
ena of high parton density QCD. Starting from an initial condition which contains free parameters 
we have numerically solved the nonlinear evolution equation, restricting ourselves to the point 
6 = 0. The parameters have been fitted to the F2 data j3], and the resulting approximate solution 
is displayed in Fig. [T] 

The 6-dependence of the solution is restored using the ansatz: 



2 The Master Equation 




(2.1) 
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Figure 1: N{b = 0) versus r 
r_L Qsix). 



where k, is related to the 6 = solution 

k{x, r^) = - ln(l - N{r^, x,b = 0)). (2.2) 
The Gaussian form for the profile function in impact parameter space was assumed i.e. 



1 6^ 
-exp{ 



TT /?2 

" ^proton 



-'^proton 



(2.3) 



where Rproton =3.1 GeV^ ^ refers to the radius of the target proton. The dipole-proton cross- 
section (from eq. ()2.4p of j4j) is given by: 



crdipoie{r±,x) = 2 d%± N{r±,b±, x) . 



(2.4) 



Another popular saturation model was proposed by Golec-Biernat and Wiistoff [2] which we will 
denote by GBW. The following dipole cross-section is assumed to be: 



^{r±,x)dipoie = o-o[l - exp{- 



4RI 



)] 



(2.5) 



with Rl{x)[GeV ^] = {^)^- The values of the parameters which were determined by fitting to 
DIS data at HERA for x"< 0.01, are: cxq = 23 mb, A = 0.29 and xq = 3 x 10"^. The r± 
dependence is taken as Gaussian, which leads to a constant cross-section ctq for large r± (or small 
Q"^) i.e. for "soft" interactions. 

The dipole cross-section of the GBW saturation model and N are closely related. Both models 
include the effects of gluon saturation, preserve unitarity, and describe the physics associated 
with "long distances". Whereas has some support from QCD, adipoie of the GBW model has 
more the character of a phenomenological model. 

Unlike the GBW (Jdipoiei the dipole cross section obtained from the solution of the BK equation is 
not saturated as a function of x. This emanates from the integration over h. With the assumed 
Gaussian profile function it leads to a logarithmic growth with decreasing x. The Froissart- 
like behaviour adipoie ln^(l/x) [131 crucially dependent on the fact that the large impact 
parameter behaviour of the dipole cross-section is exponential, rather than Gaussian. 
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Starting from our dipole cross sections we obtain our master equation for the hadron-proton 
cross-section 

(7H-proton{x) = J d^r±\^H{r±)\^ (ydipoie{r±, x) , (2.6) 

where iPh{'^±) represents the wave function of the hadron which scatters off the target proton. 
The form taken for ipH{r±) is discussed in the next section. 

For both saturation models f (j2.4j) and (j2.5|) ) the energy dependence of the hadron-proton cross 
section enters only through x-dependence of the dipole cross section, the latter being adjusted 
or constructed to describe DIS data of the F2 structure function. The way in which the the 
x-dependence of the dipole cross section determines the energy dependence of the hadron cross 
sections is strongly influenced by the 6-slope. Note, however, that there is no x in hadron-hadron 
collisions. In order to relate x to the energy of the process we will need to introduce an additional 
nonperturbative scale, denoted below as Qq. 



3 Details of Calculation 
3.1 Hadronic wave function 

There is no established method for calculating hadronic wave functions within the framework of 
QCD. The Heidelberg group Dosch et al.jllj using the stochastic vacuum model, have calculated 
hadronic cross-sections, after making an ansatz regarding the form of the hadronic wave function 
ipH{f)- We adopt their ansatz and utilize hadronic wave functions of similar shape to those used 
in 



Based on experimental evidence of the flavour dependence of hadronic cross-sections, which de- 
crease with increasing number of strange quarks, the authors of J3] hypothesized that the cross 
sections depend on the sizes of the hadrons in the process. 

For the hadron transverse wave function we take a simple Gaussian form, the square of the wave 
function is given by 



\^M{rA_)? = -^exp{-^) (3.1) 

where Sm is a parameter related to the meson size. We have used S*^ = 1.08 fm and Sk = 
0.95 fm. These Sm were found from experimental values for the electromagnetic radii, namely, 
Rt^ = 0.66 ± 0.01 fm and Rk = 0.58 ± 0.04 fm For meson wavefunction of the form ()3.1|) 



The proton's wave function squared is given by 



\Mri±,r2^)\'' = ^^exp(- "^^ +/^^ ) (3.2) 



where Sp = 1.05 fm, which corresponds to Rp = 0.862 ± 0.012 fm For proton wavefunction 



of the form ^ Sp = J^Rp. 
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For meson-proton scattering, the meson is treated as a quark-antiquark pair (i.e. a colour dipole), 
and therefore the calculation follows that of DIS, i.e. the interaction of a colour dipole with a 
proton target, with the meson wave function replacing that of the virtual photon. However, for 
the scattering of a baryon projectile, we represent the baryon as constituted of two colour dipoles, 
one dipole formed around two quarks, and the second dipole from the center of mass of these two 
quarks to the third quark in the baryon. Generally speaking the parameter Sp in ()H.2|1 can be 
different for these two dipoles. For example, in the non relativistic additive quark model (AQM) 
we expect Sp for the first dipole to be larger by the factor 4/3 compared to the second dipole. 
In ()3.2|1 . for simplicity, they are choosen to be identical. It will be shown below that the hadron 
cross sections are mostly determined by the saturation domain where the sensitivity to the wave 
function and, in particular, to the choice of Sp is quite weak. 

In AQM, the | ratio between vr-p and p-p cross-sections is due to quark counting. In our model, 
at very high energies when the dipole cross section is independent of the dipole size, the predicted 
ratio is |, given by the number of dipoles in the pion relative to those in the proton. In our 
approach the high energy interaction is blind with respect to the flavours of the interacting 
quarks, and the ratio (Jtot /(^tot = 1/2, seems to disagree with the data. Experimentally, at an 
energy of y/s = 20 GeV, crf^^/afot ~ 0.6. We expect the secondary Regge trajectories to give 
a smaller contribution to K^p interaction, as there are no resonances in the s channel of this 
reaction, the same is also true for proton-proton scattering. The predicted ratio aj^t ^/<^tot = 1/2 
is in reasonably good agreement with the experiment data. 

3.2 Method of Calculation 

All the parameters in N{r±, b±, x) were taken from the fit of 

F2{Q',x) = ^^a^*^\Q\x) (3.3) 

made to the experimental DIS data, (see |lj for details). In the DIS case the variable x is well 
defined in terms of and py^, in the hadronic case we redefine x to be x = where s denotes 
the energy squared in the center of mass system of the hadrons, and Q\ is a parameter which we 
adjust to be compatible with the data. The value of Q\ is determined by the longitudinal part 
of the wave function, for which at present we do not have a reasonable model. In general we 
would expect the scale Qo to increase with increasing hadron masses, and thus vary from hadron 
to hadron. In this study we will fit this parameter separately for each projectile hadron. 

In the colour dipole picture (which is equivalent to two gluon exchange), one can hopefully only 
reproduce the asymptotic energy dependence i.e. the Pomeron contribution. At lower energies 
where most of the data for meson baryon scattering is available, there are also contributions from 
secondary trajectories. The evaluation of these contributions is beyond the scope of our model. 
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4 Comparison of Our Model Predictions with Data 



Our model contains a parameter Qq, (which can be considered as a scale factor) from the definition 
of X = in analogy with the variable in DIS, x ~ The energy dependence of the 

hadronic cross-sections can be adjusted by choosing a suitable value for this variable. We found 
that the value Ql = 3.5 x 10^'^ GeV"^, gives good agreement with the data for crtot{K^p). 
See Fig. |21 (a). The K~^p channel was chosen as it is exotic, having no resonances in s channel 
and therefore by duality the secondary Regge trajectory contributions are small. By replacing the 
wave function of the meson in Eq. ljH.lj) by that of the vr meson, we obtain atot for the n~p. We 
display our prediction and the data in Fig. 121(a). We have chosen to show o'toti'^^p) compared to 
data, as there are more data in this channel than in the ir'^p channel. The cross-sections a^+p and 
(Ttt-p only differ due to the contribution of the secondary trajectories. The parameters used for 
atotiT^-p) are Ql = 1 x lO'^ Gel/^ and a Regge contribution of 27 x (^)-°-^^ mb (sq = I GeV^) 
which has a smaller residue than that suggested by the DL model. 

(a) (b) 



28 - 




18 - 



200 400 600 10" 10= 10= 

Plab ("SV) Vs (GeV) 

Figure 2: (a) K^p and ir'p total cross-sections. The full line is the prediction of our model and the 
dotted line using the Golec-Biernat Wiistoff dipole. 

(h) pp and pp total cross-sections. The full lines are the model predictions with a Gaussian 
profile. The dashed line is the result of using a Ki profile, and the dotted line using the 
Golec-Biernat Wiisthoff dipole. 
Data compilation from Ref. 

We note that the predicted energy dependence for crtotij^^v) is niore moderate than the data, and 
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at an energy of piab = 400 GeV, we underestimate the experimental data by approximately 7 %. 

We compare our results with those of the GBW model, by replacing adipoie in Eq- f)^-6|) . by the 
GBW dipole (j(a;, r^)rfjpoie given in Eq. ()2.5|) . We wish to stress that the GBW saturation model 
was formulated for, and applied to DIS reactions j2|. The responsibility of extending the model 
to hadronic interactions is ours. The GBW model does not contain any explicit 6-dependence, 
but one can consider the constant cross-section to be the result after integration over the impact 
parameter (with a sharp impact parameter cutoff put into the exponent in (2.5)). It is clear that 
at asymptotic energies (where one can neglect the contribution of secondary trajectories), the 
GBW cross-section for atot{Mp) < 23 mb, and atot{Bp) < 46 mb, as we have one (two) dipoles 
interacting with the proton target for Mp (Bp) scattering. In Fig. [21 (a) we show the result (dotted 
line) obtained using the GBW form for the dipole Eq. ()2.5|) for the K^p channel, with Qo ^'^^^^ = 
0.27 GeV^, and using the same form of the hadronic wave function as discussed in Section (3.1). 
For 7r~p channel we used the almost maximal possible dipole cross section, put Q"^ = b x 10~^ 
GeV^, and add a secondary trajectory with 11 x (^)~°'^^ mb. Due to the arguement presented 
above, the GBW model cannot be adjusted to these data at all. 

For dtotipp) (see Fig. El (b)) the energy dependence predicted by the model is not as steep as 
the experimental data, yielding a value for atotaiipp) ~ 65 mb (instead of 72 mb) at Tevatron 
energies, i.e. a deficit of 10 %, however, this is over a much wider energy range than in the Ti~p 
case. For atotai{pp) (where data is only available over a narrower range of energy) we achieve a 
very good reproduction of the experimental data, this is displayed in Fig. 121(b). For the pp and pp 
channels we take Ql = 0.03 GeV'^, and following |9j have a Regge contribution of 98.4 x 



-0.45 



■SQ' 



mb for atotipp) and 56 x (^) mb for atotipv)- 

To get a handle on the theoretical uncertainties of our treatment we have also calculated <Jtot{pp), 
using a profile function 

Si,) . ^Ak.(^) (4.1) 

-Tlproton -^proton -^prolan 

which corresponds to the Fourier transform of the " dipole" form factor in the momentum transfer 
representation: 

F"dipoie"it) = ^2 -. — ■ (4.2) 

The result of our calculations with the "dipole" form factor for atot{pp) is shown by the dashed 
line in Fig. 121(b). For this calculation we took Qq = 0.06 GeV^. For the comparison with the F2 
data jl], a value of R^roton = 4.46 GeV'"^ was used. 

We repeat the same procedure for the pp channel for the GBW model as we did for atot{K^p) 
explained above, now with parameters Qo ^"^^^^ = 0.7 GeV'^^ and a Regge contribution of 20 x 
(^)~°'^^ mb, (adjusted to the data), and taking the same form for the baryon wave function. The 
results are shown dotted line in Fig. |21(b). 

We also calculate the forward slope of the elastic cross- section i.e. B, which is defined as 

da da 

dt dt it^o 
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and is related to the sizes of the particles participating in the reaction. B = Bq + B' where 



B 



Jd\^\^PHir^)\ni Nir^,b^,x) db] 

(7 tot 



< 



hi > 



(4.3) 



In Fig. 01 we display B = Bq + B' with Bq = 7.8 GeV^^. Bq is related to the formfactors of 
the hadrons. Its value was chosen with an eye on the data, and is close to that used by Schuler 
and Sjostrand [T7] . 

The results we obtain are disappointing, but understandable and demonstrate the weak point of 
our model viz. the assumption of the oversimplified form for the impact parameter dependence of 
the amplitude (see Appendix). The elastic slope (unlike atot) is sensitive to the b± distribution, 
and our assumption that the major contribution comes from small values of b± is obviously wrong. 
We will expand on this difficulty in Section 5 and the Appendix. 

Using the relation 



(^elastic 



[(^tot 



16 n B 

we check our model's predictions for aeiastic for the pp channel, where the data extends to high 
energies. The model produces a very good description of the elastic pp cross-section (see Fig. 
m^a)). In aeiastic, the deficiency in the energy rise of atot is to some extend compensated by the 
inadequacy of B. 
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Figure 3: pp forward elastic slope. Data compilation from Ref. lib]/ 




Figure 4: For pp scattering; (a) elastic cross-section (h) the ratio Rd- 
Data compilation from Ref. 



We find it of interest to investigate liow close we are to a black disc picture for dipole-proton 
scattering, which we do by utilizing the Pumplin bound |19j, which follows from Unitarity con- 
siderations and can be written as: 

D CTeiastic + c^diff I d'^r _Ld%±\ilJH{r±)\'^ {r _l, b_L, x) 1 
-Kd = = < - . 

O"tot <7 tot ^ 

In Fig. Ilfb) we display the ratio Rn and the experimental data. The model's predictions agrees 
with the data, and suggests that even at Tevatron energies we are 20% away from the black disc 
limit of |. 



5 Discussion and Conclusions 

Our treatment has two parameters: 

(a) i?^ which is taken from the electromagnetic radius of the hadron (following the Heidelberg 
prescription); 

(b) the parameter Ql, (introduced after Eq. ()3.3|) . is adjusted by comparing with the data for the 
different channels. It is worth mentioning that, as expected, the obtained values for Qq reflect 
the mass hierarchy of the projectile hadrons. 

In addition, we require the contribution of secondary trajectories at lower energies. 

We would like to emphasis that our approximate solution to the BK equation ji4| is obtained for 
impact parameter b = 0, and then an ansatz is made regarding the b dependence of the profile 
function S(b). In the original fit to the F2 data |3], S(b) was taken to be a Gaussian, (which is 
equivalent to assuming that the dependence upon t, the momentum transfer squared, is exponen- 
tial). It is this Gaussian shape of the profile function which produces a cross-section with a ln(i) 
dependence. We have shown that taking a "dipole" behaviour in the t representation we would 
have a profile function S{b) ~ ^Ki{^), and the resulting cross-sections would asymptotically 
have a In^(^) dependence. A reflection of this fact is that the energy dependence of the cross 
section (see dashed line of Fig. ^ is much steeper, which is more in accord with the data. Since 
the original fit to the DIS F2 data was made with a Gaussian profile in b space, for consistency 
one should redo the fit with a different profile. This is a task for the future. We summarize our 
main results: 

• We obtain a reasonable value of the total cross-sections. 

• The predicted energy increase is too moderate. Both the value of the cross section and 
the energy dependence depend on the assumed 6-dependence of the profile function. For 
illustration we have compared with the saturation model of j2] : because of the sharp cutoff 
in b built into this model, the energy dependence is even weaker. This emphasizes the need 
to improve, in all saturation models the 6-shape of the profile function. 

• The difference between ax+p and a.^p at high energy needs an explanation, which cannot be 
answered within the framework of our model. 



10 



• The slow increase of the slope B, is the result of having an almost black disc picture for 
dipole - proton scattering [20] and not a consequence of a particular form for the 6-profile. 

• = tends to |, which is the black disc limit. We reproduce the experimental 
values for this ratio, which we consider as a success of the model. 

In general the approach works better than one would expect, even in the region of long distances. 
The results obtained are not very sensitive to the input parameters of the projectile hadron, 
namely the wavefunction and the parameter Q^. Our use of the saturation models has some 
predictive power, provided we have enough information about projectile wave function. 

Nevertheless, there are serious shortcomings in applying the dipole picture with the concept of 
saturation to hadronic cross sections. The crucial feature seems to be the 6-dependence of the 
dipole cross section which needs further investigation. On a more fundamental level, it is not 
clear at all, how to relate the dipole picture to the additive quark model: at present these two 
approaches look almost orthogonal to each other. 
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7 Appendix 

Our results can be easily checked in the asymptotic limit where we have (for simplicity of argu- 
ment the following discussion is presented at some fixed r±, while r± integrations are implicitly 
assumed) : 

aMp = 2vr y ° dbl = 2nbl{x) (7.1) 
where bl{x) is a result of our ansatz on the 6_l behaviour viz. 

(rd^pole^p{r^,x■,b^) = 2(1 - ^7 2) 

with 

n = In [1 - N{r±,x;b± 
bl{x) can be calculated from the equation 

Q{rj_,x]b_L = bo{x)) 
2 
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I.e. 

blix) = In [i ln(l - N{r^,x;b^ = 0)] 

where 0"*°* = 2'7ibl{x) and cr*p = Anb^x). The factor 2 between the cross-sections is due to the 
fact that we have two dipoles in the proton and one in the meson. 

The elastic slope is given by Eq. ()4.H|) 

, blix)/2 blix) 

^ = -JW) ^ ^ ^ ^ 

The energy dependence of B' can now be calculated 

dB' ^ 1 dbljx) ^ 1 da%\x) 
d ln(l/x) Ad ln(l/x) 16n d ln(l/x) ' ^ ' ' 

j From this formula it is clear that we have an increase of B with - which is slower than that of 
cTpp*(x), and obviously not in accord with the experimental data. Though the above argument 
was presented for the Gaussian profile, similar conclusions would be obtained for an alternative 
profile function. 
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